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Abstract 

A calculation based on flat spacetime symmetries shows how there can be two 
quantum phases. For one, extreme phase change determines a conventional classical 
trajectory and four-momentum, i.e. mass times four-velocity. The other phase occurs 
in an effective particle state, with the effective energy and momentum being the rate 
of change of the phase with respect to time and distance. A cosmic ray proton moves 
along a classical trajectory, but exists in an effective particle state with an effective 
energy that depends on the local gravitational potential. Assumptions are made so 
that a cosmic ray proton in an ultra-high energy state detected near the Earth was 
in a much less energetic state in interstellar space. A 300 EeV proton incident on the 
Earth was a 2 PeV proton in interstellar space. The model predicts such protons are 
in states with even more energy near the Sun than when near the Earth. 
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1 Introduction 

In quantum mechanics, a particle's classical trajectory maximizes or minimizes the phase 
change along the trajectory. And the way the particle state's phase changes with location 
determines the particle's momentum and energy; e.g. the momentum is proportional to the 
change in phase per unit distance. 

In Ref. [1], hereafter referred to as T, the phase that determines the trajectory and the 
phase that determines energy and momentum were found to differ in a way depending on 
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the strength of the gravitational field. Thus one of these phases sends the particle along a 
classical trajectory while the other determines the particle state's energy and momentum. It 
is shown in this article that, for sufficiently high energy particles, even a weak gravitational 
potential can produce dramatic changes to the particle state compared to the state in a null 
potential. 

In I a well-known method [2J of obtaining the quantum fields of free massive particles 
is generalized. The field is constructed as a sum over particle creation and annihilation 
operators. The method relies on the group theory of spacetime symmetries and the general- 
izations involve translations. For free particles the momentum is unchanged by translation; 
that unfaithful representation (rep) assigns unity to all translations. By generalizing to a 
nontrivial rep, momentum can change with translation, so the motion in fiat spacetime is in 
general curved. 

Having two distinct phases, 'dual phases,' derives from a well known property of trans- 
lations: any translation preserves all coordinate differences since any displacement can- 
cels upon subtraction. The same coordinate differences occur when spacetime is translated 
through one displacement 6x and the particle state is translated through a second displace- 
ment e. 

The process of constructing the quantum field constrains the displacement e of parti- 
cle states and their creation and annihilation operators. Suppose spacetime undergoes a 
Poincare transformation. Any such transformation is equivalent to a boost-rotation combi- 
nation A followed by a displacement through 6x, which gives the coordinate transformation 
X — s> Ax + 6x. 

When spacetime is transformed by {A,6x), particle states and the operators transform 
by (A, e), where the displacement e is now a free parameter. In the construction of the 
quantum field, an expression is found for e that involves a second rank tensor M, which is 
shown in I to be related to the gravitational metric tensor f/^,^, where /x, z/ G {1,2,3,4} = 
{x,y,z,t} indicate Minkowski coordinates. One must know Qf^,^, A, x, and 6x to determine 
the displacement e of the particle state. 

Perhaps unexpectedly, the particle state displacement e is not homogeneous in the space- 
time displacement 6x; there are terms proportional to x. The particle state is translated 
even when spacetime is not and the transformation applied to the particle states depends on 
location whereas the transformation applied to spacetime and the quantum field does not 
depend on location. 

This puzzle is resolved as follows. Let y denote the spacetime parameters of the particle 
states. Since the coefficient of y in the phase is the momentum, the terms in e proportional to 
y contribute to an effective momentum p which replaces the eigenmomentum p in the phase, 
p-y—>-p-y={p + de/dy) ■ y. Thus the different transformations applied to the particle 
states of a given momentum p at different locations produce a single, common plane wave 
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with an effective momentum p. The effective (four-)momentum tfierefore depends on A, and 

on Qfj^i, tlirougli M. Wc contend tliat tlie effective energy of tlic single, common plane wave 
is the energy observed by instruments that measure the energy dumped into the atmosphere 
by a cosmic ray. 

Since e depends on the Lorentz transformation A and the trajectory does not, A is free 
to vary. To see this, picture the trajectory as a sequence of translations in a given reference 
frame, i.e. 6xs not As. We choose the given frame to be determined by the distribution of 
Cosmic Microwave Background radiation. But to get to that given frame, one could start 
with any initial frame and transform to the given frame with a suitable A. Since the initial 
frame is arbitrary, and since e depends on A, one must select an initial frame to determine 
the effective momentum. 

We choose the initial frame to be the rest frame of the particle, see Fig. 1. The choice is 
based largely on its effect on the cosmic ray spectrum, but it should be noted that choosing 
the rest frame to be special is grounded in basic physical intuition. 

With the rest frame as the initial frame, A is a Lorentz transformation L[p) taking the 
rest momentum k — {0, 0, 0, m} to the trajectory momentum p. Then effective momentum 
p depends on the trajectory momentum p and the gravitational metric tensor gf^j,. The 
trajectory momentum p could be obtained by time-of flight methods and we assume that 
general relativity gives the gravitational metric tensor. It is shown in I that the trajectory 
momentum p is consistent with the gravitational metric tensor g^i,. 

Additional flexibility occurs by allowing the spacetime parameters x for the quantum 
field and the trajectory to differ from the spacetime parameters y of the particle states. 
The construction of a quantum field involves the creation and annihilation operators, not 
the particle states directly. This leaves the spacetime of the particle states y free to differ 
from X by any Lorentz transformation A. In order to avoid having the spatial momentum p 
decrease when the energy increases it is assumed that the spacetime of particle states y is 
the time inversion of the spacetime x of the quantum field and the trajectory; A is the time 
inversion transformation. 

With these assumptions, one finds that, for a null gravitational potential, = 0, the 
effective momentum p is just the trajectory momentum p — p. We are taking 'momentum' 
to mean the four-momentum, so one component is the total energy, p* = mc^ + E, where c 
is the speed of light and m is the particle mass. 

For the purpose of illustration here, consider only high energy particles, —07^ >> 1, 
where gamma ■jmc^ = mc^ + E and E is the trajectory energy obtained perhaps by a time-of 
-flight experiment. One flnds that the effective particle state energy E is approximated by 
the simple expression, 

E f» -4072^; . 
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Since 7 ~ E/{mc^), the effective energy E in the potential (j) is proportional to the cube of 
the trajectory energy E. 

Referenced to interstellar space, the potential at the Earth's surface is roughly (pE ^ 
— 1.06 X 10~^, unitless because the potential is reduced by a factor c^. The conventional 
gravitational increase in kinetic energy from interstellar space to the Earth's surface is just 
Egrav = —(pEf^c^ = 10 cV. Such Conventional contributions are negligible here and are 
ignored. For protons, the above approximation for the effective energy E is valid for 7 >> 
= 10^, i.e. for a trajectory energy E much higher than 10^^ eV, a value near the 
current maximum beam energy in accelerators. 

For example, a cosmic ray proton detected at the earth with an effective particle state 
energy of 3 x 10^° eV, the energy of the so-called 'OMG particle', [5] would have had a particle 
state energy in interstellar space of about 2 x 10^^ eV, i.e. about 5 orders lower. Since the 
potential in interstellar space vanishes, the trajectory energy in interstellar space is also 

2 X 10^5 eV. 

Cosmic rays with energies exceeding 10^° eV are very rare, so for this discussion let 

3 X 10^'' eV be the upper limit for the effective particle state energy at Earth, implying an 
upper limit of 2 x 10^^ eV for the trajectory energy which is also the upper limit for the 
effective particle state energy of a proton in interstellar space. 

A detailed discussion of the implications of this result to cosmic ray acceleration and 
transport mechanisms is beyond the scope of this article. Nevertheless, some comments 
seem apparent. The energy in interstellar space 2 x 10^^ eV is close to the limit often cited 
for acceleration of cosmic rays by the interstellar shock waves of supernovae remnants. Also, 
the energy 2 x 10^^ eV is well under the GZK cutoff of about 6 x 10^^ eV, [HE] so the proton 
would not produce pions when interacting with Cosmic Microwave Background photons. 
Finally galactic magnetic fields with observed magnitudes of 2 x 10^^° T [6l [3 E] would 
harness a 2 x 10^^ eV proton with an orbit radius of about one parsec, so the acceleration 
could take place in the Galaxy. 

The model makes specific assumptions in order to explain the energies of cosmic rays 
detected on or near the Earth's surface. Thus the cosmic ray experiments influence the 
fundamental assumptions of the model and are part of the basis for the predictions of the 
model. The main result is the dependence of particle state energy and momentum on local 
gravity. The predictions of the model include the cosmic ray spectra in interstellar space 
and at the Sun in Fig. 2. 



2 DUAL PHASES 



5 



2 Dual Phases 

In I [Ij a quantum field is constructed by modifying a well-known procedure [2J giving the 
fields of massive free particles. For convenience the construction is briefly discussed here. 

The quantum field of a particle species with nonzero mass m can be constructed as a 
linear combination of the operators that create or annihilate the single particle states of the 
particle. Let il^i{x) denote the quantum field and let a^j^jf) and aj.("p*') be the annihilation 
and creation operators that remove or add to a multiparticle state a single particle state of 
(four-)momentum p and z-component of spin a. The fourth component of the momentum, 
the energy, can be found knowing the three spatial components "p^ and the mass m. One has 

M^) = 11 d^~^ui^{x, 'f)a„{'f) + d'fvi„{x, 'f)al{'f) , (1) 

where the -us and f s are the coefficient functions. 

Clebsch- Gordon coefficients connect quantities that transform differently with the same 
rotation. The coefficient functions u and v connect the operators a and with the field ip. 
These quantities transform differently with a Poincare transformation, one (the operators) 
with unitary reps and the other (the fields) with nonunitary reps. So the -us and f s act in 
much the same way as Clebsch-Gordon coefficients. 

A single particle state for a given momentum is a plane wave that gains a phase factor 
upon translation. Translations are included in Poincare transformations of the operators 
a and and the operators are multiplied by the same phase factor as the particle states. 
But the transformation of the field ipi{x) does not introduce any phase factors. Thus, the 
coefficient functions u and v must contain compensating phase factors to adjust for the phase 
factors in the operators. 

The coefficient functions u and v keep these compensating phase factors even when no 
transformations are applied. This is why a quantum field is a sum over plane waves. 

The particle states have phases and the coefficients u and v have phases. For free fields, 
the phases of the particle states and the phases of the coefficients u and v are the same 
within a sign. For the quantum fields that respond to forces in I, the phases of the particle 
states and the phases of the coefficient functions u and v differ. 

To get fields that respond to forces in I, one allows the operators a and to transform 
with (A, e) when the field ipi{x) transforms with the Poincare transformation (A, 5x). The 
operators transform by the same Lorentz transformation A followed by a translation along a 
possibly different displacement e. The basis for allowing this is the invariance of coordinate 
differences under translations: any translation preserves all coordinate differences. And 
coordinate differences determine the spacetime intervals whose invariance is postulated in 
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special relativity. Since e can be any displacement, it can depend on A, x, and 5x, i.e. e is 
an arbitrary function e(A, x, 6x). 

To show how this goes, consider a scalar annihilation field 

^{x) = J (i^T? "jf )a("j?) . (2) 

For other spins and for details, see I. Apply the transformation (A, 6x) to the field ip and 
(A, e) to the operator a, leaving the coefficients unchanged. One has 



ij{Ax + 6x)= f (f^u{x, p^)e'^P"(^'"'^"\/^^a(Aj) , (3) 



p. 

where p* is the time component of the four-momentum. 

Since the parameter p in ([2]) is a dummy variable we can replace it with Ap. Then using 
d?Ap = d^jf{Apy/p^, see Ref. [9], and applying ([2]) with Ax + 6x in place of x, one has 

.3^ i^P)\ ~ ~ 



ip(Ax + 6x) = / d-^jt , u(Ax + 6x, Ap)a(Ap) , (4) 
J 



implying with ([3]) that 



^^^^-u(Aa; + 5x, A?) = u{x, -f)e'^P<^^-'^-) . (5) 



pt 



The coefficient u on the left is evaluated at Ax + 5x and the u on the right is evaluated at x. 
By considering a special case of ([5]) one can have the ms on both sides evaluated at x = 0. 
Consider ([5]) with A = 1 and 5x = —x. This yields 

m(x, t) = u{0, p^)e-^P-(i'-'-^-) . (6) 

This and its modification with Ax + 6x in place of x can be substituted in One gets 



^ (^P)* ^^Q^ ^ ^^^^ -^^gjAp-[e{A,x,fc)-A<:(l,a;,-x)+e(l,Ax+5x ,-Ax-5x)] _ ^^^'j 

Compare this with ([5]). Now both us are evaluated at x = 0. 

The exponential in cannot depend on x or 6x, because nothing else in does. The 
function e{A,x,6x) is thereby constrained. One can show that the displacement e(A,x, 5x) 
must have the following form, 

e'^(A, X, 6x) = e^(A) - A^M^x'^ + M'^^Alx" + M'^Jx" , (8) 
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where M = M{x) is an arbitrary second rank tensor field defined over spacetime x and M' 
= M{Ax + 6x). The next section explores some consequences of ([H]). 

We assume e'^(A) = because we are interested in the part of e that is linear in x or 6x. 
Note that one recovers e = 6x when the field M is the identity, Mj^ = b^. 

The displacement e(A, x, 5x) depends on the event x, so one must construct the field 
event-by-event, as is emphasized in I. The meaning of M can be found by considering the 
classical trajectories of a particle described by the quantum field. This is considered next. 

With expression ([8]) for e, one finds by ([6]) that 

p') = m(0, p')e*P-^^ . (9) 

Thus the coefficients u are approximately 'plane waves' in regions where p and M change 
slowly. By ([2]), is a sum over these 'plane waves'. 

Maximizing or minimizing the change in phase 5B = p ■ Mbx determines the classical 
trajectory. It is shown in I that extreme phase change occurs for Mbx proportional to p. 
Thus there is a quantity 5r such that 

^ = M,"^. (10) 

But, from the start, the momentum p is in the Wigner class with constant positive mass and 
positive energy. The magnitude of p is the mass, = — m^, and j9* > m > 0. By (ITU]) we 
can write = — as 

^"^^;^^'^^ = -i' (11) 

where rjap is the flat spacetime metric of the spacetime upon which the field is constructed. 
Define the second order tensor g^^, 

g,. = Vaf3M^M^. (12) 

Then, by combining the preceding three equations, we can write the equation = — as 

Sxf'Sx" 

for a displacement 6x along the classical trajectory. It is shown in I that g^i, can be taken 
to be the gravitational metric tensor, making r the proper time. In what follows, we assume 
that gfj^i, is known from general relativity and use f|T2|) to restrict M. Note that M is not 
uniquely determined. At least the sign of M is not determined since f|T2|) is quadratic in M. 

It is at this point in I that one invokes a nontrivial representation of translation applied 
to momentum p. With this assumption the trajectory of a particle is in general curved even 
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though the spacetime x is flat, because translating the momentum p along a displacement 
changes the momentum in accordance with the nontrivial representation of translation. In 
general relativity, by contrast, at least in its geometric interpretation, [10] a flat spacetime 
implies there are no gravitational forces and particles move in straight lines at constant speed 
in the absence of other forces. Here and in I and before that in Ref. [11], particles can move 
along curved paths in flat spacetime. 

Note that it is the momentum p that determines the trajectories by maximizing the phase 
p ■ Mx. Since observing the trajectory by measuring the time to travel a known distance, 
or some other time-of-flight method, would result in the momentum p, the momentum p is 
called the 'trajectory' momentum to distinguish it from another momentum found below in 
Sec. El 

In this section, the quantum field is shown to be a sum over coefficient functions pro- 
portional to phase factors exp {ip ■ Mx) . Extreme phase change directs the particle to follow 
classical trajectories, as is shown in I. In the next section, the expression ([8]) for the displace- 
ment e is considered, giving rise to a different phase with another momentum, an effective 
momentum that takes part in interactions. 

3 Cosmic Rays 

In this section we interpret the formula ([8]) for the displacement e of the operators and the 
particle states. The discussion leads to the momentum of particle states which is shown to 
depend on gravity. Assumptions are made so that the ultrahigh energy states of some cosmic 
rays detected in Earth's gravity would exist in much lower energy states in interstellar space. 

Let A indicate the event at xa- Quantities in an initial reference frame are indicated by 
'0'. The initial coordinates and initial momentum are xq and po- The initial coordinates of 
event A are Xq. After a transformation (A, 5x) the quantities are x -\- Sx, p, and x^ + 6x. 

We can infer some properties of particle states from the way operators transform. Single 
particle states transform in the same way as the creation and annihilation operators that add 
or remove them from multiparticle states. Thus given a transformation (A, 6x) of spacetime 
Xo, we know that the single particle state with eigenmomentum pq removed by the operator 
a(po) undergoes the same transformation (A, e^) as the operators. Here, the displacement 

indicates that e is evaluated using ([S]) at event A. 

The operator a(po) removes the same particle state wherever the quantum field is con- 
structed, but the transformation applied to the operator and therefore to the particle state 
changes from event to event. Thus the operators are not functions of xq, but the transfor- 
mations applied to the operators are functions of Xq. 

One thing that we cannot infer from the transformation of the operators is what spacetime 
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the transformation of the particle state is apphed to. Any Poincare transformation is a 
symmetry operation for any spacetime coordinate system. It may be that the spacetime yo 
for a single particle state differs from the spacetime Xq of the trajectory and quantum field. 
Let us assume that the transformed spacetimes differ by a Lorentz transformation A, 

x = Xy, (14) 

where A could be an inversion. 

Thus we deduce that a single particle state created by the operator a(po) is a plane wave 
proportional to exp {±ipo ■ yo). Particle states exp {±ipo ■ yo) and quantum fields ip{xo) = 
J d^pou{xo,jf'Q)a{jfQ) both have dependences on Jf^ and therefore on Pq. We infer that 
spacetimes Xo and yo must transform with Lorentz transformations just as po transforms, so 
that spacetime scalar products po ■ xq and po ■ yo are preserved. 

By Sec. [2l see ([3]) for the operator transformation, an initial single particle state 
exp {±ipo ■ yo) transforms with the Poincare transformation (A, e^) when spacetime under- 
goes the transformation (A, 6x). The superscript A indicates that the displacement e depends 
on event A. The displacement in general differs from the displacement e*" associated with 
a nearby event C. 

Having different transformations act at different events motivates a closer look at how 
a plane wave realizes Poincare transformations. The plane wave assigns a phase factor 
exp {±ipo ■ yo) to an event with coordinates yo. The phase factor at any event, say A, by 
itself realizes a representation of the Poincare group. The successive transformations (Ai, e^) 
followed by (A2, e^) applied to the phase factor at the event A yield 

where one recognizes the law for successive Poincare transformations, (A2, e^)(Ai, e^) = 
{A,A,,A2et + e^). 

When just one transformation is applied to a plane wave, the phase factor at each event 
undergoes the same transformation. When there is a distinct transformation (A, e^) at each 
event A, a suitable generalization is to apply the distinct transformation at each event to the 
phase factor at that event. 

Therefore, the transformations (A, e) are applied event-by-event to an initial plane wave 
exp {ipo ■ yo), yielding for the phase factor at event A, 

where the last expression follows from ([8]) and f|T4|) . p = Apo, y = Ayo, Mq is the tensor in 
the initial reference frame and M is in the transformed frame with Mq = A~^MA and M = 
AMoA-\ 
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Recall that M is related to the gravitational metric tensor g^j,^, see (|T2l) . For many 
situations with a particle moving in weak gravitational fields, the gravitational field changes 
little over a region of space confining the particle for a short time. One can treat M as 
constant over fairly large regions, large on the scale of the relevent portion of the quantum 
field. And on such a scale, the change of the momentum p is often small due to a weak 
gravitational force. Thus we can treat both M and p as constants on a scale much larger 
than the scale of the quantum field. 

Then the coefficient of y"^ in the phase of the exponential in (ITSl) does not depend on 
the event A. This means that a unique, common plane wave is formed by the process. The 
momentum of the common plane wave is an 'effective momentum' p given in 



The two momenta p and p are equal when A = 1 or M = 1. 

Keeping track of indices and displaying the flat spacetime metric r] yields an expression 
for the effective momentum. One has 



where summation over repeated indices is understood. Note that in general p is not the 
trajectory momentum p = Apo that one would observe by measuring the time for a particle 
of mass m to travel a known distance along the trajectory. 

Since the effective momentum p is the momentum of the plane wave obtained by trans- 
forming the particle state of momentum pq, we assume that this plane wave forms the effective 
particle state in the transformed frame. Since the plane wave with the effective momentum p 
is the set of phase factors at each event in spacetime that describes the state of the particle, 
we conclude that for experiments that measure the momentum and energy of cosmic rays by 
the transfer of momentum and energy to other particles, it is the effective momentum and 
energy p that is recorded. 

By making suitable assumptions for M, A, and A we can apply the expression for the 
effective momentum (fT6|) to cosmic rays so that they have less effective energy traveling in 
interstellar space (M = 1) than they have when detected in the gravitational potential of 
the Earth (M 7^ 1). 

We only consider protons as primaries. If other particles are primaries then the results 
below need to be adjusted. 

Consider a moving proton well-separated from other protons. The relevant portion of 
the proton quantum fleld can be conflned to a tube by combining particle states with eigen- 
momenta spread out over a suitable range to conform with the Heisenberg Uncertainty 
Principle. 



p.y = p.y + p. (_M + A-^MA)Xy = p{l - MX + A' 



MAX) ■ y . 



p- = / 5| - r]''%rM;X; + r]''Pr]f,,{A-'MA)lX; 



(16) 
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Spacetime outside the tube need not be translated. We assume that the external space- 
time remains in a fixed frame that we call the 'given frame'. We are most interested in 
so-called ultrahigh energy cosmic rays with energies of 10^® eV or more, moving at nearly 
the speed of light relative to the Earth and Sun. So we choose the given fiat spacetime 
reference frame to be a frame with the Earth and Sun moving with speeds negligible with 
respect to the speed of light. A suitably universal reference frame is provided by the distri- 
bution of Cosmic Microwave Background (CMB) radiation. The observed dipole anisotropy 
of the CMB implies the Solar System is moving at a speed of 370 km/s = 0.00123c [121 IE] 
with respect to the CMB reference frame. Thus the CMB reference frame is suitable for the 
purposes here. We take the given frame to be the CMB reference frame. 

Then the motion can be pictured as a succession of translations through displacements 
Sxi applied to the region of spacetime inside the tube. Since the particle is confined to the 
tube, the resulting sequence of quantum fields is just the same as if the translations were 
applied to the whole of spacetime. 

Since the motion of the proton is described as a sequence of translations with dis- 
placements in various spacetime directions, the accelerations of the proton are described 
by translations only and without the need for rotations or boosts. Here is where the nontriv- 
ial representation of translation allows and directs the proton's accelerations. No Lorentz 
transformations are applied. 

Now there is a problem: By (fT6l) . each effective momentum pi depends on an as- yet- 
unspecified Lorentz transformation Aj. Thus there is a sequence of some as-yet-unspecified 
initial reference frames that are translated to the given frame with the as-yet-unspecified 
transformations Aj. 

Lacking any reason to prefer one frame over another, we need to make an assumption 
before we can calculate the effective momenta pi. 

Assumption: Each initial frame that determines each Lorentz transformation Aj is the 
rest frame of the particle. 

Then each inital momentum is the momentum at rest, k = {0,0, 0,m}, and the 
Lorentz transformation Aj is the transformation L{pi) taking k to momentum of the particle's 
trajectory pi in the CMB reference frame; see Fig. 1. We have 



Aik = L{pi)k = Pi . 



(17) 



From here on the sequence index is dropped, e.g. Aj — A and Pi p, etc. 
There are many ways to transform k to p. For definiteness, let L be 

Ll{p) = 6l + {l + ^)-'m~Yp', 

L\ = L'l = m'^p"- and L4 = 7 = mT^p^ , 
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where m is the mass of the particle, i.e. the proton's mass. 

We turn now to M. For simphcity we assume a spherically symmetry, diagonal gravita- 
tional metric tensor g^^. For trajectories in a weak gravitational field with the gravitational 
potential 0, < and | |<< 1, one has [2 



Qf,^ = diag{g^^, g^^, g^^, gu) = diag(l - 20, 1 - 20, 1 - 20, -1 - 20) . (19) 

By f|T2|) . one choice for M is diagonal, 

M; = diag(M,, M,, M,, Mt) = diag(l - 0, 1 - 0, 1 - 0, 1 - 0, 1 + 0) , (20) 

where terms of second order in are dropped. 

The weak field gravitational potential at a point Q in space is a sum over sources, 

*=-i:^. (21) 

where G is the universal gravitational constant, is the mass of the source s, is the 
spatial distance from s to the point Q, and c is the speed of light. 

As it happens, the rotation curves for disk galaxies are fiat, so the orbital speeds of stars 
in the disk are approximately independent of their distance from the Galactic center. There 
is only a small difference between the inner and outer disk kinetic energies of circularly 
orbiting objects (per unit mass), i.e. the potential in the disk is the same as the potential 
beyond the disk. Therefore, it is reasonable to assume that the gravitational potential of the 
Galactic disk vanishes. The gravitational potential in the Galactic disk in interstellar space 
far from any star is assumed to vanish. 

Having made assumptions for M and A, it remains to consider A, which relates the 
spacetime for the trajectory and the spacetime for the effective particle states. If we assume 
the two spacetimes are the same or opposite, i.e. A = ±1 in f|T^ . then one can have the desired 
reduction in effective energy for the cosmic rays. But then the effective three-momentum 
p , the spatial part of the four-momentum p, would not behave properly. 

Heuristically, one can argue as follows. The lowest order term M = 1 cancels out of 
( IT6|l . The next order term has —0 for and -|-0 for Mt, see ( 120|) . making this part of M 
proportional to the matrix diag{ — 1, — 1,— 1,+1}. It must be that MX can be proportional 
to 1 (the delta function = diag{+l, +1, +1, +1}), since the first term in parentheses in (fT6l) 
is the delta function. A spatial inversion or a time inversion for A combined with the part 
of M linear in yields a matrix MX proportional to 1, which can be expected to behave 
reasonably. One suspects that A could be an inversion. 
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One finds that assuming A to be a time inversion produces a spatial momentum that 
increases when the energy increases, as one would expect. Thus, we assume that A is the 
time inversion 

A^ = diag(+1,+1,+1,-1). (22) 

With the above assumptions for M, A, and A one can obtain a suitable expression for the 
effective momentum p. 

By (|T6l) . (|T7l) . (1201) . and (l22i) we have the effective momentum 

p'^=p^(l-47V), (23) 

where only the lowest order terms in (p are kept, indicates the spatial part of the effective 
four-momentum, k G {1,2,3}, and = = mc^ + ^ is the effective total energy of the 
particle state, including the rest energy. Recall that p is the trajectory four-momentum 
observable by time-of-fiight measurements. 

Ultra-high energy cosmic rays are observed with more energy than expected. By fl25]) the 
effective energy is more than the energy of the trajectory since < 0. Because experiments 
determine the energy of the particle state, the measured cosmic ray energy is the quantity 
E = p^ — mc^ with the total energy given in Eq. fl23l) with evaluated at the Earth. For 
the accuracy needed here, the thickness of the atmosphere and the altitude of Earth orbiting 
satellites can be neglected. By (|2T1) . we find the gravitational potential (p to be —1.06 x 10~^ 
at the Earth's surface due to the Sun and the Earth. 

Experiments that measure the energy of cosmic rays have succeeded in pushing the 
observed spectrum to 'ultra-high' energies E > 10^^ eV. Some of the data [15] - [21] is 
collected in the spectrum of Fig. 2 labeled 'Earth'. Not all available data is included 
because the sketch is already busy and the discussion here involves the overall properties of 
the spectrum upon which all data agree. 

The expected energy spectrum of the protons in interstellar space (I.S.) can be obtained 
from the energies of cosmic rays observed by Earth-based experiments. One applies the chain 
rule to obtain the predicted scaled flux in interstellar space using (l23l) . 

/ rfA^ \ _ f dN \ dlnE 
\ dlnE J I.S. \ dlnE J Earth dlnE ' 

where {dN/ dlnE) Earth is the scaled flux observed by Earth-bound experiments, E = p^ — mc^ 
and E = p^—mc^. The resulting energy spectrum is plotted in Fig. 2 and labelled 'Interstellar 
Space'. 

Just as fl2^ and deduce the energy spectrum in interstellar space from the observed 
earth-based spectrum, one can reverse the process and predict the cosmic ray spectrum 
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for any weak potential (f). The gravitational potential at the surface of the Sun is (j) — 
—2.12 X 10~^, with (f) — in interstellar space as before. The spectrum is plotted in Fig. 2 
and labelled 'Sun'. 

Note that the product of the energy E and the flux dN/dE is plotted in Fig. 2, EdN/dE 
— dN/dlnE. We call this the 'scaled flux'. The scaled flux is closely related to the number 
of particles. Since the same number of particles (per unit area per unit time per steradian) 
are involved in each spectrum, this fact explains the alignment of the three spectra. 

One sees from Fig. 2 that the spectra overlap up to about 10^^ eV, where the spectra 
diverge and at the highest energies the spectra differ by many orders of magnitude. The 
predicted interstellar spectrum may be compared with the astrophysics of accelerating and 
transportating the protons. Such a comparison lies beyond the scope of this paper. It may 
be possible to detect cosmic rays striking the Sun where the highest energy cosmic rays in 
Fig. 2 should have energies more than a hundred times larger than those striking the Earth. 

Experimental confirmation of the predicted cosmic ray spectra in interstellar space and 
at the Sun's surface would provide support for the assumptions and explanations presented 
in this article. 
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A Problems 

1. A proton approaches the Earth radially with a trajectory energy of E = 2 x 10^^ eV as 
measured in the Cosmic Microwave Background reference frame. Find the Lorentz contracted 
diameter of the Earth in the rest frame of the proton. 

2. Show that if M satisfies Eq. ( 1121) . then AM also satisfies the equation with the same g^u. 
Here A is an arbitrary Lorentz transformation, i.e. a sequence of rotations and boosts, but 
no translations. 

3. Show that the expression for e in ([H]) makes the phase factor in ([7]) independent of x and 
6x. 

4. The spectra in Fig. 2 straighten out at high energy, (i) Fit each spectrum at high energy 
to a power law, 

where a and n are constants. Note that Fig. 2 plots EdN/ dE. (ii) Considering an energy Eq 
in the power law portion of each spectrum, use the values of a and n from (i) to estimate the 
number of cosmic rays with energy Eq or higher and the intensity passing through one square 
meter in one second from each incident steradian. (iii) For Earth, compare the intensity with 
the Solar Constant which is about 1400 W/m^. (iv) Use the results of (ii) to explain why 
the high energy spectra for the Earth and the Sun align as drawn in Fig. 2. 

5. Using the approximate effective momentum (|23|) . find dN/dE as a function of dN/dE, E, 
and (p. 

6. The OMG particle deposited 3 x 10^° eV of energy into the Earth's atmosphere. [3] How 
much energy would the particle have deposited into the Sun's photosphere? Assume the 
OMG particle is a proton. 

7. Get M and E for a general diagonal metric (^^j^. Assume that M is diagonal and reduces 
to fl20|) in weak fields. Then use the Schwarzchild metric in isotropic coordinates to get the 
energy of an incident proton at the event horizon of a collapsed object. What is the value 
of the energy at the event horizon for an incident proton that had an energy of 2 x 10^^ eV 
in interstellar space far from the collapsed object? 

8. Suppose SN A and SN B are supernovae with shock waves. Assume a proton is accelerated 
from shock wave A and collides with a molecule in the atmosphere of a planet B ahead of 
shock wave B. A tiny fraction of the copious energy deposited into the atmosphere of planet 
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B is used to eject a proton from planet B into interstellar space with a very low velocity. 
The proton is accelerated by shock wave B and is then incident on a planet A ahead of shock 
wave A. The proton collides with a molecule in the atmosphere of planet A, releasing a lot 
of energy. Some of that energy releases protons into interstellar space to be accelerated by 
the shock wave A and so on. Obtain an expression for the energy gained in one such cycle. 

9. Find the ratio of a component of the particle state momentum p^^ to the trajectory 
momentum component as a function of (= p* = m'j) and the weak gravitational 
potential (p using Eq. (I23l) . Assuming the particle is a proton, graph the energy ratio p'^/p^ 
as a function of logp^ for from 1 GeV to 10 TeV on the Earth's surface. Could experiments 
with currently operating proton accelerators distinguish p^ and p^7 

10. Find the effective mass rh of the effective particle state with momentum p in (l23l) as a 
function of both the gravitational potential (p and the trajectory gamma 7 = p^/mc^. Plot 
m vs 7 in interstellar space, at the Earth's surface, and on the surface of the Sun. 

11. Derive an expression like fl23l) except assume that the spacetime for particle states is the 
negative of the trajectory spacetime, i.e. A = —1. What happens to the spectra in Fig. 2? 
What happens to the effective three-momentum p ? What happens to the effective mass rh 
in Problem 10? 
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Figure 1: The Lorentz Transformation A. Top: In the rest frame of the proton, the Earth is 
Lorentz contracted. Bottom: In the frame at rest with respect to the distribution of Cosmic 
Microwave Background photons, the Earth is nearly at rest and the proton has a trajectory 
four-momentum p, as measured by time-of-flight and mass. However the effective particle 
state of the proton has a different four-momentum p. The trajectory and the particle state 
arise from plane waves with different phases. It is p, not p, that is deposited into the Earth's 
atmosphere when the particle state changes upon interaction with the atmosphere. 
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Figure 2: Cosmic Ray Spectra, Dependence on Local Gravity. The spectrum labeled 'Earth' 
plots the experimentally determined scaled flux incident on Earth's atmosphere. [I5]- [21] The 
flux is scaled by the kinetic energy E, i.e. dN/dlnE = EdN/dE. The spectrum 'Interstellar 
Space' is the expected spectrum where the gravitational potential vanishes. For the most 
energetic cosmic rays, the predicted Interstellar Space spectrum is remarkably less energetic 
than the spectrum observed on Earth. At the Sun the magnitude of the gravitational poten- 
tial is greater than at the Earth and the predicted cosmic ray energies at the Sun are higher 
than those detected at the Earth. 



